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ABSTRACT 

We propose a methodology to perform a self-consistent analysis of the physical 
properties of the emitting gas of HII galaxies adequate to the data that can be ob- 
tained with the XXI century technology. This methodology requires the production 
and calibration of empirical relations between the different line temperatures that 
should superseed currently used ones based on very simple, and poorly tested, photo- 
ionization model sequences. 

As a first step to reach these goals we have obtained simultaneous blue to far red 
longslit spectra with the William Herschel Telescope (WHT) of three compact Hll 
galaxies selected from the Sloan Digital Sky Survey (SDSS) Data Release 2 (DR2) 
spectral catalog using the INAOE Virtual Observatory superserver. Our spectra cover 
the range from 3200 to 10500 A, including the Balmer jump, the [On] AA 3727,29 A 
lines, the [Sill] AA 9069,9532 A doublet as well as various weak auroral lines such as 
[Oin] A 4363 A and [Sin] A 6312 A. 

For the three objects we have measured at least four line temperatures: T([0iii]), 
T([Siii]), T([Oil]) and T([Sli]) and the Balmer continuum temperature T(Bac). These 
measurements and a careful and realistic treatment of the observational errors yield 
total oxygen abundances with accuracies between 5 and 9%. These accuracies are 
expected to improve as better calibrations based on more precise measurements, both 
on electron temperatures and densities, are produced. 

We have compared our obtained spectra with those downloaded from the SDSS 
DR3 finding a satisfactory agreement. The analysis of these spectra yields values of line 
temperatures and elemental ionic and total abundances which are in general agreement 
with those derived from the WHT spectra, although for most quantities they can only 
be taken as estimates since, due to the lack of direct measurements of the required 
lines, theoretical models had to be used whose uncertainties are impossible to quantify. 

The ionization structure found for the observed objects from the + /0 2+ and 
S + /S 2+ ratios points to high values of the ionizing radiation, as traced by the values 
of the "softness parameter" r\ which is less than one for the three objects. The use 
of line temperatures derived from T([0iii]) based on current photoionization models 
yield for the two highest excitation objects, much higher values of rj which would imply 
lower ionizing temperatures. This is however inconsistent with the ionization structure 
as probed by the measured emission line intensities. 

Finally, we have measured the T(Bac) for the three observed objects and derived 
temperature fluctuations. Only for one of the objects, the temperature fluctuation is 
significant and could lead to higher oxygen abundances by about 0.20 dex. 

Key words: galaxies: fundamental parameters - galaxies: starburst - galaxies: abun- 
dances - galaxies: temperature - ISM: abundances - Hll regions: abundances 
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1 INTRODUCTION 

Hn galaxies are low mass irregular galaxies 
with, at least . a recent episode of vi olent star 
formation jMelnick. Terlevich fc EggletorJ Il985l : 

iMelnick. TerleviclT^^MolesI ^8^1 concentrated in a 
few parsecs close to their cores. The ionizing fluxes orig- 
inated by these young massive stars dominate the light 
of this subclass of Blue Compact Dwarf galaxies (BCDs) 
which show emission line spectra very si milar to those of 
giant extragalactic H ll regions (GEHRs; ISargent fc Searlel 
Il970fc iFrenchl Il980l) . Therefore, by applying the same 
measurement techniques as for Hll regions, we can derive 
the temperatures, densities and chemical composition of 
the interstellar gas in this type of generally metal-deficient 
galaxies. In some cases, it is possible to detect in these 
objects, intermediate-to-old stellar populations which have 
a more uniform spatial distribution than the bright and 
youn g stellar populations associate d with the ionizing clus- 
ters fechulte-Ladbeck et alj |l998fl . This older population 
produces a characteristic spectrum with absorption features 
which mainly af fect the hydrogen recombination emission 
lines (lDiazHl988T) . that is the Balmer and Paschen series in 
the spectral range of interest. In some cases, the underlying 
stellar absorptions can severely affect the ratios of Hi 
line pairs and hence the determination of the reddening 
constant (C(H/3)). They must therefore be measured with 
special care (see discussion in f£JJ. 

A considerable number of the blue objects observed 
at intermediate redshifts seem to have properties (mass, 
R e , velocity width of the em i ssion l ines) similar to local 
Hu g alaxies llKoo et al.l 11994 Il995l iGuzman et"aH Il996l . 
119981). In particular, those with CT<65kms _1 follow the 
same a — Lb and Lhb relation as seen in Hu galaxies 
iMelnick et alJl20qd: iMas- Hesse et al.ll2003l : pTerlevich et"ai] 
120031: ISieeel et al.ll2005l) ~ Similar conclusion is drawn from 
recent studies on Lyman Break galaxies that also suggest 
that strong narrow emission line galaxies might have been 
very common m the past (e.g. IPettini et all l200l 1200 ll : 
lEllison et all 1200 ll) . To detect possible evolutionary effects 
like systematic differences in their chemical composition, ac- 
curate and reliable methods for abundance determination 
are needed. 

This is usually done by combining photoionization 
model results and observed emission line intensity ratios. 
There are several major problems with this approach that 
limit the confidence of present results. Among them: the 
effect of temperature struct ure in multiple-zone models 
iPerez-Montero fc Dfazl 12003:); the presence of tem pera- 
ture fluctuations across the nebula JPeimbertlEooil : col- 
lisional and density effects on ion temperatures (Luridi- 
ana et al. 1999, Perez- Montero & Diaz 2003); the pres- 
ence of neutral zones affecting the calculation of ion- 
ization correction factors (ICFs; Peimbert et al. 2002); 
the ionization s tructure not adequate l y rep roduced by 
current models iPerez-Montero fc Dfazl 120031) ; ionization 
vs. matter bounded zones, affecting the low ionization 
lines formed in the outer parts o f the ionized regions 
JCastellanos. Diaz fc Terlevich! Eo53) . On the other hand, 
the understanding of the age and evolutionary state of Hll 
galaxies require the use of self-consistent models for the 
ionizing stars and the ionized gas. However, model com- 
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puted evolutionar y sequences show impor tant differences 
with observations iStasiriska fc Izotovl |200 Ij) , including: (a) 
Hell is too strong in a substantial number of objects as com- 
pared to model predictions; (b) [Om]/H/3 vs. [Oll]/H/3 and 
[Om]/H/3 vs. [Ol]/H/3 are not well reproduced by evolution- 
ary model sequences in the sense that predicted collision- 
ally excited lines are too weak compared to observations; (c) 
there is a large spread in the [Nil] /[On] values (more than 
an order of magnitude) for galaxies with the same value of 
([Oll] + [Om])/H/3 in the metallicity range from 8 to 8.4 (see 
e.g. Perez-Montero & Diaz 2005). 

Substantial progress toward solving the problems listed 
above has to come from the accurate measurement of weak 
emission lines which will allow to derive [On], [Sn] and 
[Sill] temperatures and densities allowing to constrain the 
ionization structure as well as Balmer and Paschen dis- 
continuities which will provide crucial information about 
the actual values of temperature fluctuations. It is pos- 
sible that these fluctuations produce the observed differ- 
ences between the abundances relative to hydrogen de- 
rived from recombination lines (RLs) and collisionally ex- 
cited lines (CE Ls) when a constant electron tempera- 
ture is assumed l)Peimbertlll967l : IPeimbert fc Costerolll969t 
IPeimberdlTimh . These discrepancies have been observed 
in a good sample of objects, such as galactic Hn re- 
gions (e.g. lEsteban et alJ 120041 : lOarci'a-Roias et al.1 1200.4 
120061 and reference s therein), Hn regions in the Mag- 
ellanic Clouds (e.g . IPeimbert et~aH 120001 IPeimbertl 120051 
iTsamis et alJ 120031. and references therein), extragalactic 
Hn regions (e.g. IPeimb ert fc Peimbertl2003t IPeimbert et alJ 
120051: iGuseva et al.ll200 J and referenc e s therein) and plan- 
etary nebulae ( e.g.TRubin et al.ll2002l: IWesson et al.l 20051: 
iLiu et alJlioOrj lLTull2006l : IPeimbert fc Peimbertl 12009. and 
references therein). Likewise, there are relatively recent the- 
oretical works that study the possible causes of these dis- 
crepancies in abundance determinations us ing photoioniza- 
tion models of different complexity (e.g. |s tasiriska| 120051 : 
Ijamet et all2005l : ITsamis fc Peouignoll2005l . and references 
therein). 

Unfortunately most of the available starburst and Hn 
galaxy spectra have only a restricted wavelength range (usu- 
ally from about 3600 to 7000 A) and do not have the ade- 
quate S/N to accurately measure the intensities of the weak 
diagnostic emission lines. Even the Sloan Digital Sky Survey 
(SDSS) spectra (Stoughton et al. 2002) do not cover simul- 
taneously the [On] AA 3727,29 and the [Sm]A9069A lines. 
We have therefore undertaken a project with the aim of ob- 
taining a database of top quality line ratios for a sample 
that includes the best objects for the task. The data is col- 
lected using exclusively two arm spectrographs in order to 
guarantee both high quality spectrophotometry in the whole 
spectral range from 3500 to 10500 A and good spectral and 
spatial resolution. In this way we will be in a position to 
vastly improve constraints on the photoionization models 
including the mapping of the ionization structure and the 
measurement of temperature fluctuations about which very 
little is known. 

In this first work we present observations of three Hn 
galaxies selected from the SDSS. Details regarding the se- 
lection of the objects as well as the observations are given 
in Section 2. Section 3 presents the results including line 
measuring techniques. The methodology for the derivation 



of gaseous physical conditions and elemental abundances is 
presented in Sections 4 and 5 respectively. The discussion of 
our results, including a detailed comparison with the SDSS 
data is presented in Section 6. Finally, Section 7 summarizes 
the main conclusions of our work. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Object selection 

SDSS constitutes a great base for statistical studies of the 
properties of galaxies. At this moment, the Fourth Data Re- 
lease 1 (DR4), the last one up to now, contains five band 
photometric data for about 18xl0 7 objects and more than 
6.7xl0 5 spectra of galaxies, quasars and stars (Adelman- 
McCarthy et al. 2006). The spectroscopic data have a reso- 
lution (R) of 1800-2100 covering a spectral range from 3800 
to 9200 A, with a single 3arcsec diameter aperture. The 
SDSS data were reduced and flux-calibrated using automatic 
pipelines. However, when we selected our objects on July 
2004, the DR2 2 was just available. Both data releases, DR2 
and DR4, contain the same type of objects observed in the 
same five photometric bands and using the same spectro- 
scopic configuration. DR2 has photometric data for over 88 
million unique objects and about 3.7x10 spectra (Abaza- 
jian et al. 2004) . All the objects belonging to one data release 
are also included in the next ones, but some objects could 
have been re-calibrated or re-observed somehow. Using the 
implementation of the SDSS database in the INAOE Virtual 
Observatory superserver 3 , we have selected from the whole 
SDSS DR2 the brightest nearby narrow emission line galax- 
ies with very strong lines. Our selection parameters were: 

Equivalent width of Ha > 50 A, 

1.2<a(Ha)<7A, 

redshift , z < 0.2 and 

Ha flux, F(Ha) > 4 x 10~ 14 erg cm- 2 s" 1 A" 1 . 

This preliminary list was then processed using BPT 
(Baldwin, Phillips & Terlevich 1981) diagnostic diagrams 
to remove AGN-like objects. 

The final list consists of about 200 bonafide bright Hn 
galaxies. They show spectral properties indicating a wide 
range of gaseous abundances and ages of the underlying stel- 
lar populations (lL6pezll2005ri . 

From this list, the final set was selected by further re- 
stricting the sample to the largest Ha flux and highest S/N 
objects. 

From the 10 brightest of the final set we selected three 
HII galaxies to be observed in our one single night observ- 
ing run. This final selection was made based on the relative 
positions of the sources in the sky allowing to optimize the 
observing time. The journal of observations is given in ta- 
bleland the photometric characteristics of the objects are 
listed in tabled Subsequently we have used the explore tool 4 
implemented in the DR3 5 (Abazajian et al. 2005), which was 



1 http://www.sdss.org/dr4/ 

2 http://www.sdss.org/dr2/ 

3 http://ov.inaoep.mx/ 

4 http:/ /cas. sdss.org/astro/en/tools/cxplorc/ 

5 http://www.sdss.org/dr3/ 
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Table 1. Journal of observations. Right ascension, declination and redshift were obtained from the SDSS. 



Object ID 


spSpcc SDSS 


RA 


Dec 


z 






Exposure (s) 


SDSS J002101. 03+005248.1 


spSpec-51900-0390-445 


5.254297 


0.880038 


0.098 


1 


X 


1200+2 X 2400 


SDSS J003218.60+150014.2 


spSpec-51817-0418-302 


8.077479 


15.003949 


0.018 


1 


X 


1200+2 X 2400 


SDSS J162410.11-002202.5 


spSpec-52000-0364-187 


246.042122 


-0.367378 


0.031 


1 


X 


1200 + 3 x 1800 



Table 2. SDSS photometric magnitudes obtained using the DR4 explore tools. 



Object ID 


u 


g 


r 


i 


z 


SDSS J002101. 03+005248.1 


17.56 


17.35 


17.51 


16.98 


17.45 


SDSS J003218. 60+150014.2 


17.04 


16.49 


16.53 


16.74 


16.65 


SDSS J162410.11-002202.5 


17.07 


16.46 


16.91 


16.80 


16.74 



Table 3. WHT instrumental configuration 



Spectral range 


Disp. 


FWHM 


Spatial res. 


(A) 


(Apx- 1 ) 


(A) 


("px- 1 ) 


blue 3200-5700 


0.86 


2.5 


0.2 


red 5500-10550 


1.64 


4.8 


0.2 



available at the time of analysis, to extract again the three 
object SDSS spectra for comparison purposes. 

2.2 Observational details 

The blue and red spectra were obtained simultaneously us- 
ing the ISIS double beam spectrograph mounted on the 4.2m 
William Herschel Telescope (WHT) of the Isaac Newton 
Group (ING) at the Roque de los Muchachos Observatory, 
on the Spanish island of La Palma. They were acquired on 
July the 18th 2004 during one single night observing run 
and under photometric conditions. EEV12 and Marconi2 de- 
tectors were attached to the blue and red arms of the spec- 
trograph, respectively. The R300B grating was used in the 
blue covering the wavelength range 3200-5700 A (centered at 
A C = 4450A), giving a spectral dispersion of 0.86 A pixel 1 . 
On the red arm, the R158R grating was mounted provid- 
ing a spectral range from 5500 to 10550 A (A c = 8025 A) and 
a spectral dispersion of 1.64 A pixel -1 . In order to reduce 
the readout noise of our images we have taken the observa- 
tions with the 'SLOW CCD speed. The pixel size for this 
set-up configuration is 0.2 arcsec for both spectral ranges. 
The slit width was ~0.5 arcsec, which, combined with the 
spectral dispersions, yielded spectral resolutions of about 
2.5 and 4.8 A FWHM in the blue and red arms respectively. 
All observations were made at paralactic angle to avoid ef- 
fects of differential refraction in the UV. The instrumental 
configuration, summarized in table [3] was planned in order 
to cover the whole spectrum from 3200 to 10550 A provid- 
ing at the same time a moderate spectral resolution. This 
guarantees the simultaneous measurement of the Balmer 
discontinuity and the nebular lines of [On] AA 3727,29 and 
[Sill] AA 9069,9532 A at both ends of the spectrum, in the 
very same region of the galaxy. A good signal-to-noise ratio 
was also required to allow the detection and measurement of 
weak lines such as [Om] A 4363, [Sn] AA4068, 6717 and 6731, 
and [Sin] A 6312. 

Several bias and sky flat field frames were taken at the 
beginning and at the end of each night in both arms. In addi- 
tion, two lamp flat fields and one calibration lamp exposure 
were performed at each telescope position. The calibration 



lamp used was CuNe+CuAr. The images were processed 
and analyzed with IRAF routines in the usual manner. The 
procedure includes the removal of cosmic rays, bias substrac- 
tion, division by a normalized flat field and wavelength cali- 
bration. Typical wavelength fits were performed using 30-35 
lines in the blue and 20-25 lines in the red and polynomials 
of second to third order. These fits have been done at 117 
different locations along the slit in each arm (beam size of 10 
pixels) obtaining rms residuals between ~0.1 and ~0.2pix. 

In the last step, the spectra were corrected for atmo- 
spheric extinction and flux calibrated. For the blue spectra, 
four standard star observations were used, allowing a good 
spectrophotometric calibration with an estimated accuracy 
of about 5%. Unfortunately, only one standard star could be 
used for the calibration of the red spectra. Nevertheless, af- 
ter flux calibration, in the overlapping region of the spectra 
taken with each arm, the agreement in the average contin- 
uum level was good. 



3 RESULTS 

The WHT spectra of the observed galaxies with some of 
the relevant identified emission lines are shown in Figs. Q |2] 
and OH Each fi gure is split into two panels, showing the blue 
(upper panel) and red (lower panel) spectral ranges. 

The emission line fluxes were measured using the 
SPLOT task of IRAF and are listed for the three observed 
galaxies in tables 2] El and [(J Column 1 of each table shows 
the wave length and the name of th e measured lines, as re- 
ferred in iGarcfa-Roias et al] (l2004fl . The observed emission 
line fluxes, F(X) (in units of H/3 flux = 1000) with their cor- 
responding errors, are presented in column 2. The measured 
equivalent widths (EW) are listed in column 3. 

We have used two different ways to integrate the 
flux of a given line: (1) in the case of an isolated line 
or two blended and unresolved lines the intensity was 
calculated integrating between two points given by the 
position of the local continuum placed by eye; (2) if 
two lines are blended, but they can be resolved, we 
have used a multiple gaussian fit procedure to estimate 
individual fluxes. Following IConzalez-Deleado et al] 

6 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is op- 
erated by the Association of Universities for Research in Astron- 
omy, In. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 
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Figure 1. WHT spectra of SDSS J002101. 03+005248.1 in the rest frame for the two observed ranges. 
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Figure 2. WHT spectra of SDSS J003218.60+150014.2 in the rest frame for the two observed ranges. 
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Figure 3. WHT spectra of SDSS J162410.11-002202.5 in the rest frame for the two observed ranges. 



Jl994f). I Castellanos. Diaz fc Terlevichl i2002t) and 
|p^rez^toiter^^^Dlaal20o3) ~the statistical errors associ- 



ated with the observed emission fluxes have been calculated 
using the expression oi = <r c N 1/2 [l + EW/(NA)] 1/2 ; where 
o i is the error in the observed line flux, a c represents the 
standard deviation in a box near the measured emission 
line and stands for the error in the continuum placement, 
N is the number of pixels used in the measurement of the 
line flux, EW is the line equivalent width, and A is the 
wavelength dispersion in angstroms per pixel. There are 
several lines affected by bad pixels, internal reflections 
or charge transfer in the CCD, telluric emission lines or 
atmospheric absorption lines. These cause the errors to 
increase, and, in some cases, they are impossible to quantify. 
In these cases we do not include these lines in the tables 
nore in our calculations. The only exception is the emission 
line [Sm]A9069 for SDSS J162410.11-002202.5 which is 
affected by the strong narrow water- vapor lines present in 
the A 9300 - 9500 wavelength region foiaz. Pagel fc Wilson] 
Il985l) . We have listed the value of the measurement of this 
line in table |SJ but all the physical parameters depending 
on its intensity were calculated using the theoretical ratio 
between this line a nd [Sm]A9532, 1(9069) w 2.44 x 1(9532) 
JOsterbrocldll989ri . In some cases there is an observable 
line (e.g., [Clm] AA 5517,5537, several carbon recombination 
lines, Balmer or Paschen lines) for which it is impossible 
to give a precise measurement. This might be due to a 
low signal to noise between the line and the surrounding 
continuum. This is also the case for the Paschen jump 
that could not be measured even though it was observed, 



because it was very difficult to locate the continuum at 
both sides of the discontinuity with an acceptable accuracy. 

The spectrum of SDSS J002101. 03+005248.1 presents 
very wide lines (FWHM « 7.5 A for A w 6600 A) for the ex- 
pected velocity dispersion in a low mass galaxy of this type. 
This could be due to an intrinsic velocity dispersion in this 
object, the interaction with another unobservable object or 
a projection effect on the line of sight. There are Hll galaxies 
that in fact are multiple systems, with two or more compo- 
nents, despite their "a priory" assum ption of compactness 
JZwickvlll966l ISargent fc Searlelll970l) . In some cases, these 
systems s how some evidence of interaction a mong their com- 
ponents (iTelles. Melnick fc Terlevichl 11997?) . For instance, 
IIZw40, which was fir s t clas sified as a compact emission 
line galaxy bv ISargentl lll970l) . when observed with enough 
spatial resol ution showed to be the merge of tw o separate 
subsystems fealdwin. Spinrad fc Terlevicblll982T) . As a con- 
sequence, lines that should be resolved are blended in the 
spectrum of SDSS J002101.03+005248.1. Such is the case of 
Ha and [Nil] A 6548 A emission lines. We have resorted to 
the theoretical ratio, 1(6584) ~ 3 x 1(6548), to decontami- 
nate the observed flux of Ha by the emission of [Nil] A 6548 
and to derive the electron temperature of [Nil]. 

A conspicuous underlying stellar population is easily ap- 
preciable by the presence of absorption features that depress 
the Balmer and Paschen emission lines. The upper panel 
of Fig. [I] shows an example of this effect for the Balmer 
lines (H13 to H5) on an enlargement of the spectrum of 
SDSS J003218.60+150014.2, the object that presents the 
most relevant and appreciable absorption lines. The pseudo- 
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Figure 4. Upper panel: Section of the spectrum of SDSS 
J003218.60+150014.2 taken with the WHT. The spectrum is in 
the rest frame and its spectral range is from 3725 to 4125 A. 
We can appreciate the presence of absorption features originated 
in an oldish stellar population which mainly affect the Balmer 
emission lines. Lower panel: Section of the same spectrum with a 
spectral range from 3785 to 3885 A. We have superposed the fit 
to H9 made using the ngaussfit task from IRAF (dashed-dotted 
lines). For both panels: The dashed line traces the continuum and 
the dotted lines show the pseudo-continuum used to measure the 
Balmer emission lines. 



traction of the absorption component is almost inappreciable 
for the stronger lines. In the other two galaxies, the absorp- 
tion wings in the Balmer lines are not prominent enough as 
to provide sensible results by the multi-gaussian component 
fitting. Therefore, we have doubled the error derived using 
the the expression for the statistical errors associated with 
the observed emission fluxes, (crj). 

The absorption features of the underlying stellar 
population may also affect the helium emission lines 
to some extent. However, these absorption lines are 
narrower than those of h ydrogen (see, for example, 
iGonzalez-Deleado et al.| |2005f). and therefore we cannot see 
their wings at both sides of the emission and we cannot de- 
fine a pseudo-continuum to measure the line fluxes. 

The reddening coefficient (C(H/3)) has been calculated 
assum ing the galactic extinction law of iMiller fc Mathews! 
lll972l) with R v =3.2. C(Hf3) was obtained in each case 
by performing a least square fit to the ratio between 
F(X) and F(H0) to the theoretical values computed by 
IStorev fc Hummeil lll995l) using an iterative method to es- 
timate T e and n e in each case. We have taken n e equal to 
n([Sll]). Due to the large error introduced by the presence 
of the underlying stellar population (see discussion above) 
only the four strongest Balmer emission lines (Ha, H/3, H7 
and H<5) have been taken into account. The values obtained 
for C(Hf3) and their corresponding errors, considered to be 
the uncertainties of the least square fittings, are listed in 
tables and HJ for each of the observed objects. 

The emission line intensities corrected for reddening, 
relative to H/? (1(A)), and their corresponding errors are 
listed in column 4 of tables [I] [5] and |S| The errors were 
obtained propagating in quadratures the observational er- 
rors in the emission line fluxes and the reddening constant 
uncertainties. We have not taken into account errors in 
the theoretical intensities since they are much lower than 
the observational ones. Finally, the values listed in Col- 
umn 5 of the tables indicate the fractional error in the 
line intensities calculated as explained above. These errors 
vary from a few percent for the more intense nebular emis- 
sion lines (e.g. [Oni] AA 4959,5007, [S11] AA 6717,6731 or the 
strongest Balmer emission lines) to 10-16 % for the weakest 
lines that have less contrast with the continuum noise (e.g. 
HelAA 3820,7281, [Anv]A4740 or OlA8446). For the auro- 
ral lines, the fractional errors are between ~3 and ~9 %. 



continuum used to measure the line fluxes is also shown. We 
can clearly see the wings of the absorption lines implying 
that, even though we have used a pseudo-continuum, there 
is an absorbed fraction of the emitted flux that we are not 
able to measure with an acceptable accuracy (see discussion 
in Diaz 1988). This fraction is not the same for all lines, nore 
the ratios between the absorbed fractions and the emissions 
are the same. In order to quantify the effect of the underly- 
ing absorption on the measured emission line intensities, we 
have performed a multi-gaussian fit to the absorption and 
emission components seen in this galaxy. The fitting can be 
seen in the lower panel of Fig. [I] The difference between 
the measurements of the absorption subtracted lines and 
the ones obtained with the use of the pseudo-continuum is, 
for all Balmer lines, within the observational errors and, in 
fact, the additional fractional error introduced by the sub- 



4 PHYSICAL CONDITIONS OF THE GAS 

4.1 Electron densities and temperatures from 
forbidden lines 

The physical conditions of the ionized gas, including elec- 
tron temperatures and electron density, have been derived 
from the emission line data using the same procedures as in 
Perez-Montero & Diaz (2003), based on the five- level sta- 
tistical equilibrium atom approximation in the task TEM- 
DEN, of the software package IRAF (De Robertis, Dufour 
& Hunt, 1987; Shaw & Dufour, 1995). The atomic coeffi- 
cients used here are the same as in Perez-Montero & Diaz 
(2003; see table 4 of that work), except in the case of + for 
which we have used the transition probabilities from Zeip- 
pen (1982) and the collision strengths from Pradham (1976), 
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Table 4. Relative observed and reddening corrected line intensities [F(H(3)=I(H(3)=W00] for SDSS J002101. 03+005248.1 



A (A) 



SDSS J002101. 03+005248.1 - spSpcc-51900-0390-445 



WHT 





F(X) 


-EW(A) 


/(A) 


Error (%) 


F(A) 


-EW(A) 


7(A) 


Error (%) 


3697 H17 


4+1 


0.3 


6+1 


11.7 










3704 H16+Hei 


11+1 


0.7 


15+1 


9.9 


— 


— 


— 


— 


3712 H15 


8+1 


0.5 


12+2 


15.5 


— 


— 


— 


— 


3727 [On] 6 


1190+9 


71.1 


1634+19 


1.2 


1598+15 


95.0 


1789+20 


1.1 


3750 H12 


21+2 


1.3 


29+3 


9.4 


27+3 


1.8 


30+4 


12.4 


3770 Hll 


28+2 


1.7 


38+3 


7.8 


30+3 


2.0 


33+4 


11.3 


3798 H10 


31+2 


2.0 


42+3 


6.7 


35+3 


2.3 


39+4 


9.5 


3835 H9 


49+5 


3.2 


65+6 


9.7 


57+6 


4.0 


63+7 


11.1 


3868 [Nein] 


294+6 


17.7 


388+8 


2.1 


346+6 


21.6 


382+7 


1.9 


3889 Hei+H8 


155+9 


9.7 


203+11 


5.6 


169+9 


11.2 


186+10 


5.5 


3968 [Nem]+H7 


226+10 


15.2 


290+13 


4.4 


248+14 


17.2 


271+15 


5.5 


4026 [Nii]+Hci 


12+1 


0.8 


15+2 


10.9 


10+1 


0.6 


10+1 


11.9 


4068 [Sn] 


11+1 


0.8 


14+1 


8.6 


12+1 


0.8 


12+1 


9.7 


4102 H<5 


213+6 


15.5 


265+7 


2.8 


230+7 


17.1 


249+8 


3.2 


4340 H7 


423+7 


33.6 


499+9 


1.7 


442+9 


36.5 


469+9 


2.0 


4363 [Om] 


47+2 


3.8 


56+3 


4.9 


40+2 


3.3 


42+2 


5.2 


4471 Hei 


37+3 


3.0 


42+3 


7.9 


36+2 


3.1 


38+3 


6.9 


4658 [Fein] 


10+1 


0.8 


10+1 


7.1 


11+1 


0.9 


11+1 


12.3 


4686 Hen 


6+1 


0.5 


7+1 


8.1 


4+1 


0.3 


4+1 


7.0 


4713 [Ariv]+Hci 


6+1 


0.6 


7+1 


5.6 


5+1 


0.5 


5+1 


10.5 


4740 [ArlV] 


5+1 


0.1 


5+1 


12.4 


— 


— 


— 


— 


4861 H/9 


1000+9 


97.0 


1000+9 


0.9 


1000+9 


108.1 


1000+9 


0.9 


4881 [Fein] 


— 


— 


— 


— 


4+1 


0.4 


4+1 


11.2 


4921 Hei 


6+1 


0.6 


6+1 


8.6 


8+1 


0.9 


8+1 


6.9 


4959 [Om] 


1575+13 


151.5 


1532+13 


0.8 


1518+10 


160.0 


1503+10 


0.6 


4986 [Fem] c 


9+1 


0.9 


9+1 


10.2 


10+1 


1.1 


10+1 


10.5 


5007 [Om] 


4514+26 


439.4 


4334+26 


0.6 


4583+26 


489.4 


4516+26 


0.6 


5199 [Ni] 


13+1 


1.5 


12+1 


11.2 


12+1 


1.4 


11+1 


8.1 


5270 [Fcm] a 


— 


— 


— 


— 


9+1 


1.1 


9+1 


9.1 


5755 [Nil] 


7+1 


0.9 


5+1 


7.1 


7+1 


1.0 


6+1 


6.9 


5876 Hei 


161+8 


22.3 


127+7 


5.2 


137+7 


22.0 


125+6 


5.0 


6300 [Oi] 


54+2 


8.2 


39+1 


3.3 


47+3 


8.4 


42+3 


6.1 


6312 [Sm] 


17+1 


2.6 


12+1 


4.8 


16+1 


2.8 


14+1 


6.8 


6364 [Oi] 


15+1 


2.4 


11+1 


5.2 


14+1 


2.6 


12+1 


6.8 


6548 [Nil] 


— 


— 


— 


— 


— 


— 


— 


— 


6563 Ha 


4160+39 


658.3 


2886+40 


1.4 


3236+8 


607.1 


2841+22 


0.8 


6584 [Nil] 


377+12 


59.6 


260+9 


3.3 


268+6 


50.8 


235+5 


2.2 


6678 Hei 


47+2 


7.7 


32+1 


4.6 


38+2 


7.4 


33+2 


6.0 


6717 [Sn] 


201+7 


33.8 


136+5 


3.4 


180+5 


35.4 


157+5 


3.0 


6731 [Sn] 


159+7 


27.2 


107+5 


4.3 


136+5 


26.9 


118+4 


3.6 


7065 Hei 


39+2 


7.2 


25+1 


4.8 


30+2 


6.5 


25+1 


5.8 


7136 [Arm] 


97+4 


18.5 


62+2 


3.9 


74+3 


17.5 


63+2 


4.0 


7155 [Fen] 










3+1 


0.8 


3+1 


12.5 


7254 Oi 










3+1 


0.8 


3+1 


13.1 


7281 Hci a 


10+1 


1.9 


6+1 


7.0 


8+1 


2.1 


7+1 


5.6 


7319 [On] d 


31+1 


6.2 


20+1 


3.2 


24+1 


5.8 


21+1 


2.5 


7330 [On] e 


21+1 


1.2 


13+1 


4.5 


18+1 


1.2 


15+1 


3.6 


7378 [Nin] 


6+1 


1.2 


4+1 


8.4 


4+1 


0.9 


3+1 


8.9 


7412 [Nin] 


7+1 


1.4 


4+1 


7.9 










7751 [Arm] 


23+2 


5.1 


14+1 


6.9 


22+2 


6.1 


19+1 


7.4 


8446 Oi 


28+2 


7.2 


16+1 


7.6 










9069 [Sm] 


230+11 


71.8 


119+6 


5.0 










9532 [Sin] 


476+16 


504.4 


239+9 


3.8 










I(H/3)(ergseg _1 cm -2 ) 
C(H/3) 


2.47 x 10~ 14 
0.51+0.01 


3.38 x 10" 14 
0.18+0.01 



SDSS 



"possibly blend with an unknown line; 6 [On] AA 3726 + 3729; c [Fern] AA 4986 + 4987; d [On] AA 7318 + 7320; e [On] AA 7330 + 7331. 
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Table 5. Relative observed and reddening corrected emission line intensities [F(H(3)=I(H(3)=W00] for SDSS J003218.60+150014.2. 



A (A) 



SDSS J003218.60+150014.2 - spSpec-51817-0418-302 



WHT 





F(X) 


-EW(A) 


/(A) 


Error (%) 


F(A) 


-EW(A) 


7(A) 


Error (%) 


3687 H19 


10+1 


0.6 


10+1 


8.5 










3697 H17 


13+1 


0.9 


14+1 


8.6 


— 


— 


— 


— 


3704 H16+Hei 


19+2 


1.2 


20±2 


8.4 


— 


— 


— 


— 


3712 H15 


15±1 


0.9 


16+1 


7.5 


— 


— 


— 


— 


3727 [On] 6 


1466±7 


80.0 


1573±14 


0.9 


— 


— 


— 


— 


3734 H13 


18±1 


1.2 


20±1 


7.5 


— 


— 


— 


— 


3750 H12 


20+2 


1.3 


21 ±2 


7.7 


— 


— 


— 


— 


3770 Hll 


34+2 


2.3 


37±2 


4.6 


45±3 


3.0 


51±4 


7.3 


3798 H10 


42+2 


2.8 


45±3 


5.7 


52±4 


3.6 


59±4 


6.8 


3820 Hei 


7±1 


0.1 


7±1 


14.3 


— 


— 


— 


— 


3835 H9 


61+5 


4.4 


65±5 


8.3 


59±5 


3.7 


67±5 


7.7 


3868 [Nein] 


350±9 


19.9 


372+10 


2.6 


370±11 


19.7 


418±12 


2.9 


3889 Hei+H8 


184±12 


12.5 


196+12 


6.3 


181±11 


10.7 


204±13 


6.2 


3968 [Nem]+H7 


243±12 


17.0 


257±13 


5.0 


244±15 


17.4 


272±17 


6.1 


4026 [Nn]+Hel 


Mil 


0.8 


14+1 


7.1 


15±2 


0.9 


17±2 


12.1 


4068 [Sn] 


18±1 


1.1 


19+1 


5.5 


13±1 


0.8 


15±1 


9.1 


4102 H<5 


236±9 


17.2 


247±9 


3.7 


237±8 


17.5 


261±9 


3.3 


4340 H7 


482±10 


33.8 


500±10 


2.0 


440±8 


32.5 


473±9 


1.9 


4363 [Om] 


60+3 


4.2 


62+3 


4.8 


60±4 


4.2 


65±4 


6.4 


4471 Hei 


36+2 


2.7 


37±2 


5.6 


34±3 


2.1 


36±3 


8.5 


4658 [Fein] 


9±1 


0.7 


9±1 


9.1 


10±1 


0.7 


10±1 


8.1 


4686 Hen 


13+1 


1.0 


13+1 


6.3 


Mil 


1.0 


14±1 


8.6 


4713 [Ariv]+Hei 


11+1 


0.9 


11+1 


5.5 


11±1 


0.8 


11±1 


11.4 


4740 [Anv] 


5±1 


0.4 


5±1 


12.9 


3±1 


0.3 


3±1 


11.9 


4861 Hf3 


1000±8 


89.6 


1000±8 


0.8 


1000±9 


86.7 


1000±9 


0.9 


4921 Hei 


9±1 


0.8 


9±1 


6.5 


11±1 


0.9 


llzbl 


9.6 


4959 [Om] 


1592±13 


140.7 


1582±13 


0.8 


1651±11 


138.6 


1631±11 


0.6 


4986 [Fem] c 


16+1 


1.4 


15+1 


9.1 


18±2 


1.5 


18±2 


9.2 


5007 [Om] 


4649±24 


416.4 


4607±24 


0.5 


4863±16 


412.8 


4778±16 


0.3 


5015 Hei 


23+2 


2.0 


22±2 


8.8 


24±2 


2.1 


24±2 


7.8 


5199 [Ni] 


9±1 


0.9 


9±1 


5.4 


10±1 


0.9 


9±1 


7.7 


5270 [Fcm] a 


— 


— 


— 


— 


7±1 


0.7 


7±1 


8.6 


5755 [Nil] 


— 


— 


— 


— 


— 


— 


— 


— 


5876 Hei 


131+6 


16.0 


124±6 


4.5 


125±6 


15.0 


113±5 


4.4 


6300 [Oi] 


30+1 


4.1 


28±1 


3.9 


30±1 


3.9 


26±1 


2.5 


6312 [Sm] 


25+1 


3.3 


23±1 


3.7 


20±1 


2.6 


17±1 


3.2 


6364 [Oi] 


11+1 


1.5 


10±1 


5.9 


10±1 


1.3 


9±1 


3.0 


6548 [Nil] 


41+3 


5.8 


37±2 


6.5 


37±2 


5.2 


31±2 


5.9 


6563 Ha 


3089+19 


435.5 


2848±30 


1.0 


3309±16 


462.8 


2825±27 


0.9 


6584 [Nil] 


121+5 


17.3 


112±5 


4.4 


108±3 


15.1 


92±3 


3.3 


6678 Hei 


35+2 


5.2 


33±2 


5.4 


34±1 


4.9 


29±1 


2.4 


6717 [Sn] 


188+5 


30.3 


173±5 


2.9 


185±4 


26.9 


156±3 


2.1 


6731 [Sn] 


136+3 


21.9 


125±3 


2.5 


135±3 


19.7 


114±2 


2.1 


7065 Hei 


24+2 


4.2 


22±2 


9.6 


24±2 


3.9 


20±1 


6.6 


7136 [Arm] 


102+3 


17.7 


92±3 


3.5 


91±4 


15.0 


75±4 


4.9 


7155 [Fen] 


5±1 


0.9 


4±1 


13.6 










7281 Hel a 


6±1 


1.1 


6±1 


12.0 


5±1 


0.9 


4±1 


6.3 


7319 [On] d 


27+1 


5.0 


25±1 


4.0 


24±1 


4.1 


20±1 


2.9 


7330 [Oil] 6 


21+1 


3.9 


19±1 


4.6 


18±1 


3.1 


15±1 


3.0 


7751 [Arm] 


26+3 


5.1 


23±2 


10.1 


21±1 


3.9 


17±1 


6.1 


8446 Oi 


7±1 


1.6 


6±1 


15.7 


6±1 


1.5 


5±1 


11.0 


8468 P17 


3±1 


0.8 


3±1 


12.5 










8546 P15 










7±1 


1.8 


5±1 


10.8 


8599 P14 


14±1 


3.7 


12±1 


10.7 


10±1 


2.8 


8±1 


8.6 


8865 Pll 


17+2 


5.1 


15±2 


11.2 


16±2 


1.6 


12±1 


9.4 


9014 P10 


25+3 


8.1 


21±2 


10.7 










9069 [Sm] 


251+10 


76.4 


217±10 


4.4 










9532 [Sm] 


518±43 


175.5 


445±38 


8.5 










9547 P8 


75+6 


43.6 


64±5 


8.5 










I(H/3)(crg seg -1 cm -2 ) 
C(H/3) 


1.23 x 10~ 14 
0.11+0.01 


3.17 x 10" 14 
0.22±0.01 



SDSS 



a possibly blend with an unknown line; b [On] AA 3726 + 3729; c [Fern] AA 4986 + 4987; d [On] AA 7318 + 7320; e [On] AA 7330 + 7331. 
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Table 6. Relative observed and reddening corrected emission line intensities [F(H(3)=I(H(3)=W00] for SDSS J162410.11-002202.5. 



A (A) 



3687 
3692 
3697 
3704 
3712 
3727 
3750 
3770 
3798 
3820 
3835 
3868 
3889 
3968 
4026 
4068 
4102 
4340 
4363 
4471 
4658 
4686 
4713 
4740 
4861 
4881 
4921 
4959 
4986 
5007 
5015 
5159 
5199 
5270 
5755 
5876 
6300 
6312 
6364 
6548 
6563 
6584 
6678 
6717 
6731 
7065 
7136 
7254 
7281 
7319 
7330 
7751 
8392 
8413 
8438 
8446 
8468 
8503 
8546 
8599 



H19 
H18 
H17 

H16+Hei 

H15 

[On] b 

H12 

Hll 

H10 

Hei 

H9 

[Nem] 

Hci+H8 

[Nem]+H7 

[Nn]+Hei 

[Sn] 

ns 

Uy 

[Om] 

Hei 

[Fein] 

Hen 

[Anv]+Hei 

[Anv] 

H/3 

[Fein] 

Hei 

[Om] 

[Feoi] c 

[Om] 

Hei 

[Feii] a 

[Ni] 

[Fem] a 

[Nil] 

Hei 

[Oi] 

[Shi] 

[Oi] 

[Nil] 

Ha 

[Nil] 

Hei 

[Sn] 

[Sn] 

Hei 

[Arm] 

Oi 

Hei a 

[On] d 

[On] e 

[Arm] 

P20 

P19 

P18 

Oi 

P17 

P16 

P15 

P14 



SDSS J162410.11-002202.5 - spSpec-52000-0364-187 




WHT 






SDSS 




F(A) 


-EW(A) 


/(A) 


Error (%) 


F(X) 


-EW(A) 


/(A) 


Error (%) 


4zbl 


0.3 


5±1 


10.9 










5zbl 


0.4 


6±1 


9.7 












8zbl 


0.8 


11±1 


9.2 












16±1 


1.5 


22±2 


7.6 












16±1 


1.5 


22±2 


7.8 












1066±13 


91.0 


1471zb22 


1.5 


1160±13 


96.81 


1480zbl9 


1.3 


20±1 


1.9 


28±2 


6.7 


18±2 


1.57 


23zb 2 


9.3 


24±2 


2.2 


33±2 


6.7 


24zb2 


2.02 


30zb 3 


8.3 


35±2 


3.2 


48±3 


6.7 


36±3 


2.93 


45zb 4 


8.4 


8zbl 


0.7 


llzbl 


10.9 


8±1 


0.61 


lOzb 1 


10.1 


51±4 


4.6 


69±5 


7.7 


49±3 


4.17 


61zb 4 


6.6 


322±8 


25.8 


428±11 


2.6 


359±8 


27.28 


444zbl0 


2.3 


111±9 


12.8 


186±11 


6.1 


145±9 


12.61 


178zbll 


6.3 


249±14 


21.1 


321±18 


5.7 


238zbl3 


20.46 


289zbl5 


5.3 


13±1 


1.0 


16±1 


5.9 


13±1 


1.11 


16zb 1 


6.0 


12±1 


0.9 


15±1 


8.2 


llzbl 


0.96 


14zb 1 


10.8 


207±6 


17.0 


258±8 


3.0 


212zb6 


19.74 


251zb 7 


2.9 


386±7 


34.6 


457±8 


1.8 


407zb6 


39.78 


463zb 7 


1.5 


59±2 


5.0 


70±2 


3.4 


53zbl 


5.02 


60zb 1 


2.3 


37±1 


3.1 


42±2 


3.8 


35zbl 


3.41 


39zb 1 


3.2 


9zbl 


0.8 


9±1 


8.2 


lOzbl 


0.96 


lOzb 1 


10.0 


9zbl 


0.8 


lOzbl 


7.0 


8zbl 


0.75 


8zb 1 


6.7 


11±1 


2.0 


11±1 


8.4 


lOzbl 


1.01 


lOzb 1 


7.8 


5zbl 


0.4 


5±1 


8.4 


5zbl 


0.48 


5zb 


10.5 


lOOOzbll 


100.6 


lOOOzbll 


1.1 


1000zb9 


117.31 


lOOOzb 9 


0.9 


3zbl 


0.3 


3±1 


10.3 


9zbl 


0.31 


9zb 1 


10.6 


10±1 


1.0 


10±1 


4.3 


8zbl 


0.85 


8zb 1 


8.7 


1993zbl4 


191.7 


1938zbl4 


0.7 


2005zbl2 


220.10 


1963zbl2 


0.6 


13±1 


1.3 


13zbl 


6.6 


llzbl 


1.18 


lOzb 1 


6.3 


5882±28 


565.3 


5642±28 


0.5 


6042zbl8 


659.56 


5855zbl8 


0.3 


29±3 


2.8 


28±3 


11.4 


19zb2 


2.07 


18zb 2 


10.6 


4zbl 


0.4 


4±1 


12.6 












8zbl 


0.9 


7±1 


8.0 


8zbl 


0.85 


7zb 1 


8.2 


5zbl 


0.6 


5±1 


10.2 


4zbl 


0.47 


4zb 


7.9 


3zbl 


0.4 


3±1 


8.4 












172±4 


20.4 


135±3 


2.4 


137zb3 


21.54 


113zb 3 


2.4 


42±1 


5.5 


30zbl 


2.2 


36zbl 


6.41 


28zb 1 


2.2 


25±1 


3.3 


18±1 


2.7 


19zbl 


3.31 


15zb 


3.1 


15±1 


2.0 


11±1 


3.0 


13zbl 


2.24 


lOzb 


3.3 


49±2 


6.8 


34±2 


5.1 


38zb2 


6.82 


28zb 1 


4.1 


4053zb33 


566.3 


2795±36 


1.3 


3737zbl7 


680.37 


2820zb24 


0.9 


135±3 


18.8 


93±2 


2.6 


106zb2 


19.27 


80zb 1 


1.8 


51±3 


7.7 


35±2 


5.9 


45zb3 


8.51 


34± 2 


5.9 


203zb5 


29.2 


136±3 


2.5 


169zb3 


32.21 


125zb 3 


2.0 


147±4 


21.3 


99±3 


2.7 


123zb2 


23.61 


91zb 2 


2.1 


41±3 


6.2 


26±2 


6.9 


34zb2 


7.25 


25zb 2 


6.6 


127zb4 


20.5 


80±3 


3.3 


102zb3 


22.14 


72zb 2 


2.8 


2=bl 


0.4 


1±1 


9.2 










9zbl 


1.5 


6±1 


7.3 


7zbl 


1.62 


5zb 


6.6 


35±1 


5.5 


21zbl 


3.1 


29zbl 


6.30 


20zb 1 


2.5 


28±1 


1.5 


18±1 


3.3 


24zbl 


5.16 


16zb 


2.7 


34±2 


5.8 


20±1 


6.1 


25zbl 


5.81 


16zb 1 


5.0 


9zbl 


2.0 


5±1 


7.6 


7zbl 


2.16 


4zb 


8.1 


9zbl 


2.1 


5±1 


11.0 










8zbl 


1.7 


4±1 


10.6 


5zbl 


1.67 


3zb 


8.9 


11±1 


2.3 


6±1 


7.7 


lOzbl 


3.06 


6zb 


4.9 


8±1 


1.9 


5±1 


9.7 


6zbl 


1.96 


4zb 


9.1 


20±2 


5.2 


llzbl 


10.8 










13±1 


3.5 


7±1 


11.1 


9zbl 


3.46 


6zb 1 


9.2 


12zbl 


2.7 


7±1 


9.3 


13zbl 


4.07 


8zb 1 


9.1 
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Table 6 - continued Relative observed and reddening corrected emission line intensities \F{H0)- 
002202.5. 



J(H/3)=1000] for SDSS J162410.ll- 



A (A) 




SDSS J002101. 03+005248.1 


- spSpcc-51900-0390-445 






WHT 








SDSS 






F(X) 


-EW(A) 


/(A) 


Error (%) 


F(X) 


-EW(A) 


/(A) 


Error (%) 


8665 P13 


15±1 


3.6 


8±1 


8.3 


16±1 


5.73 


10± 1 


7.5 


8751 P12 


23+2 


6.0 


12±1 


9.2 


18±2 


5.94 


11 + 1 


8.7 


8865 Pll 


33+3 


9.0 


17±2 


10.4 


24±2 


8.24 


15± 1 


10.1 


9014 P10 


34+6 


8.3 


17±3 


17.6 










9069 [Sm]' 


137+9 


66.8 


70±5 


7.0 










9532 [Sin] 


809±36 


220.0 


401±19 


4.8 










9547 P8 


91±3 


25.8 


45+2 


4.1 










I(H/3)(ergseg _1 cm -2 ) 




5.03 x 10~ 14 






8.10 x 10~ 14 




C(H/3) 




0.52+0.01 






0.39+0.01 





a possibly blend with an unknown line; b [On] 
f affected by atmospheric absorption bands. 



AA 3726 + 3729; c [Fein] AA 4986 + 4987; d [On] AA 7318 + 7320; e [On] AA 7330 + 7331; 



Table 7. Emission-line ratios used to derive electron densities and 
temperatures. 



Diagnostic 


Lines 


n([Sll]) 


I(6717A)/I(673lA) 


T([Oni]) 


(I(4959A)+I(5007l))/I(4363A) 


T([Oll]) 


I(3727A)/(I(7319A)+I(7330A)) 


T([Sm]) 


(I(9069A)+I(9532A))/I(6312A) 


T([Sn]) 


(I(6717A)+I(673lA))/(I(4068A)+I(4074A)) 


T([Nn]) 


(I(6548A)+I(6584A))/I(5755A) 



which offer more reliable nebular diagnostics results for this 
species (Wang et al., 2004). We have taken as sources of 
error the uncertainties associated with the measurement of 
the emission-line fluxes and the reddening correction and we 
have propagated them through our calculations. Electron 
densities have been derived from the [Sll] AA 6717 / 6731 A 
line ratio, which is representative of the low-excitation zone 
of the ionized gas. In all cases they were found to be lower 
than 200 particles per cubic centimeter, well below the crit- 
ical density for collisional deexcitation. We have tried to de- 
rive the electron densities from the [Anv] AA 4713 / 4740 A 
line ratio, by decontaminating the first one from the Hel 
contribution, but the derived density values had unaccept- 
able errors due to their large and sensitive dependencies on 
the emission line intensities and the errors of the observed 
fluxes. 

Several electron temperatures for each spectrum have 
been measured: T([0n]), T([Om]), T([Sn]), T([Sm]) and 
T([Nn]). The [Nil] line at 5755 A was not detected in the 
spectrum of SDSS J003218.60+150014.2 and therefore it was 
not possible to measure T([Nll]) for this object. Both the 
[On] AA 7319,7330 A and the [Nil] A 5755 A lines can have a 
contribution by direct recombination which increases with 
temperature. Using the calculated [Oiii] electron tempera- 
tures, we have estimated these contributions to be less than 
4% in all cases and therefore we have not corrected for this 
effect, although we have considered them when estimating 
the errors. The emission-line ratios used to calculate each 
temperature are summarized in table |7| 

The derived electron densities and temperatures for the 




3500 3600 3700 3800 3900 400C 

Wavelength ( A) 

Figure 5. Enlargement of the spectrum of SDSS 
J003218.60+150014.2 taken with the WHT. Its spectral 
range is from 3430 to 4000 A, and it is in the rest frame. The 
solid lines trace the continuum to both sides of the Balmer jump 
and the dashed line depict the value of its measurement. 



three observed objects are listed in columns 2, 4 and 6 of 
table |H| along with their corresponding errors. 



4.2 Balmer temperature 

The Balmer temperature depends on the value of the Balmer 
jump (BJ) in emission. To measure this value we have ad- 
justed the continuum at both sides of the discontinuity 
(A/3 =3646 A). Fig. shows an example of the procedure 
for SDSS J003218. 59+150014. 2. The contribution of the un- 
derlying population (see Section 3) affect, among others, the 
hydrogen emission lines near the Balmer jump. The increase 
of the number of lines toward shorter wavelengths produces 
blends which tend to depress the continuum level to the right 
of the discontinuity and precludes the application of multi- 
gaussian component fittings. We have taken special care in 
the definition of this continuum by using a long baseline and 
spectral windows free of absorption lines. The uncertainties 
due to the presence of the underlying stellar population (dif- 
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Table 8. Electron densities and temperatures for the observed galaxies using SDSS and WHT spectroscopy. Temperatures marked with 
asterisks have been deduced using equations from grids of photoionization models. 





SDSS J002101. 03+005248.1 


SDSS J003218.60+150014.2 


SDSS J162410.11-002202.5 




WHT 


SDSS 


WHT 


SDSS 


WHT 


SDSS 


n([Sn]) 


120±68 


77±58 


54: 


56: 


58: 


66: 


TQOiii]) 


1.25±0.02 


1.13±0.02 


1.28±0.02 


1.28±0.03 


1.24±0.01 


1.16±0.01 


TQOii]) 


1.03±0.02 


1.05±0.02 


1.35±0.04 




1.31±0.03 


1.23±0.03 


(T([Oll])>* 








1.36 






TQSiii]) 


1.31±0.05 




1.36±0.07 




1.26±0.04 




(T([Siii])>* 




1.02 




1.21 




1.06 


TQSii]) 


0.86±0.06 


0.76±0.06 


1.03±0.05 


0.92±0.06 


1.04±0.07 


1.02±0.09 


T([Nn]) 


1.19±0.05 


1.36±0.06 






1.42±0.08 




(T([Nn])>* 






1.35 


1.36 




1.23 


T(Bac) 


1.24±0.27 




0.96±0.16 




1.23±0.22 




T(H/3) 


1.24±0.31 




1.02±0.20 




1.24±0.26 




T(Heii) 


1.24±0.29 




0.98±0.19 




1.24±0.23 





densities in cm 3 and temperatures in 10 4 °K 



ferent possible continuum placements) have been included in 
the errors of the measurements of the discontinuities. They 
are actually smaller than the error introduced by the fitting 
of stellar templates. Once the Balmer jump is measured, 
the Balmer continuum temperature (T(Bac)) is determined 
from the ratio of the Balmer jump flux to t he flux of the Hl l 
Balmer emission line using equation (3) in lLiu et al.l l)200ll) : 

T(Bac) = 368 x (1 + 0.259t/ + + 3.409y 2+ ) ( J K 

where y + and y 2+ are the ionic abundances of singly and 
doubly ionized helium, He + /H + and He 2+ /H + (see O . re- 
spectively, and BJ is in ergs cm -2 s _1 A . The derived val- 
ues for T(Bac) are also listed in table |H1 



5 CHEMICAL ABUNDANCES 

Ionic and total abundances of He, O, S, N, Ne, Ar and Fe 
are obtained as detailed below and given in tables [5] and 
I1UI We have derived the ionic chemical abundances of the 
different species using the stronger available emission lines 
detected in the analyzed spectra. The total abundances have 
been derived by taking into account, when required, the un- 
seen ionization stages of each element, resorting to the most 
widely accepted ionization correction factors (ICF) for each 
species: 

± = ICF(X+*) ■ ±_ 



5.1 Helium abundance 

We have measured emission fluxes for 10 lines of Hel and 
1 of Hell, although four of them are blended with another 
emission line and two are so weak that they cannot be used 
to derive the helium abundance with the necessary accu- 
racy. We have therefore used the Hei A A 4471, 5876, 6678 and 
7065 A, and Hell A 4686 A lines to estimate the abundances 
of helium once and twice ionized respectively. These lines 
arise mainly from pure recombination, however they could 
have some contribution from collisional excitation as well 



as be affected by self-absorp tion and, if present, by under- 
lying stellar absorption (see lOlive fc Skillmanl 12001) . 1200-3. 
for a complete treatment of these effects). We have taken 
the electron temperature of [Oni] as representative of the 
zone where the He emission arises (i. e. T(Hell) ~ T([Offl])) 
and we have used the equations given by Olive & Skillman 
to derive the He + /H + value, using the theoretical emissivi- 
ties scaled to H/3 from lSmitsl lll996l) an d the expressions for 
the co llisional correction factors from iKingdon fc Ferlandl 
lll995l) . We have not taken into account, however, the cor- 
rections for fluorescence (three of the used helium lines have 
a negligible dependence with optical depth effects and the 
observed objects have low densities) and the underlying stel- 
lar population. The three galaxies show in their spectra the 
signature of the presence of WR stars by the blue 'bump' 
around A 4600 A therefore we have to take care when mea- 
suring the emission line flux o f Hen A 4686 A (see Fig. [(jjl. 
We have used equation (9) from lKunth &i Sargentl lll98St) to 
calculate the abundance of twice ionized helium. Then, for 
the total abundance we have directly added the two ionic 
abundances. 

He _ He+ + He 2+ 
IT ~ H+ 

The results obtained for each line and their corresponding 
errors are presented in tabled] along with the adopted value 
for He+/H+. 

5.2 Ionic and total chemical abundances from 
forbidden lines 

The oxygen ionic abundance ratios, + /H + and 2+ /H + , 
have been derived from the [On] AA 3727,29 A and [Oin] 
AA4959, 5007 A lines respectively. The simultaneous deter- 
mination of T([Oll]) and T([Om]) allows a more confident 
estimate of the total abundance of oxygen for which we have 
used the approximation: 

0_ _ Q+ + Q 2+ 
H ~ H+ 

Regarding sulphur, we have derived S + abundances us- 
ing T([Sn]) values and the fluxes of the [Sn] emission lines 
at AA 6717, 6731 A. In the same way, the abundances of S 2+ 
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SDSS J002101. 03+005248.1 


SDSS J003218.60+150014.2 


SDSS J162410.11-002202.5 




A (A) 


WHT 


SDSS 


WHT 


SDSS 


WHT 


SDSS 


He+/H+ 


4471 


0.086+0.006 


0.077±0.006 


0.076±0.004 


0.074±0.006 


0.086±0.003 


0.079±0.002 




5876 


0.096±0.005 


0.094±0.006 


0.096+0.004 


0.088±0.003 


0.104±0.002 


0.086+0.002 




6678 


0.087±0.004 


0.089±0.006 


0.089+0.004 


0.079±0.001 


0.094±0.005 


0.091±0.005 




7065 


0.094±0.005 


0.102±0.006 


0.086+0.009 


0.076±0.006 


0.104±0.008 


0.099+0.007 




adopted 


0.090±0.005 


0.090±0.009 


0.087+0.007 


0.080±0.005 


0.098±0.007 


0.085+0.009 


He 2 +/H+ 


4686 


0.0006±0.0001 


0.0004±0.0001 


0.0011+0.0001 


0.0013+0.0001 


0.0008±0.0001 


0.0007±0.0001 


(He/H) 




0.091+0.005 


0.090±0.009 


0.088±0.007 


0.081±0.005 


0.098±0.007 


0.085±0.009 



Table 10. Ionic and total chemical abundances, along with their corresponding ionization correction factors, derived from CELs. 
Logarithm of the S23 parameter. 





SDSS J002101. 03+005248.1 


SDSS J003218.60±150014.2 


SDSS J162410.11-002202.5 




WHT 


SDSS 


WHT 


SDSS 


WHT 


SDSS 


12+log(0+/H+) 
12+log(0 2 +/fl'+) 


7.73+0.06 


7.73+0.05 


7.15±0.06 


7.10 


7.18±0.06 


7.30±0.05 


7.86+0.02 


8.01+0.02 


7.86±0.02 


7.87±0.05 


7.98±0.02 


8.08±0.01 


12+log(0/H) 


8.10+0.04 


8.19+0.03 


7.93±0.03 


7.93 


8.05±0.02 


8.14±0.02 


12+log(S+/H+) 


5.93+0.11 


6.13+0.12 


5.80±0.06 


5.88±0.08 


5.69±0.08 


5.67±0.10 


12+log(S 2 +/H+) 


5.91+0.05 


6.55 


6.16±0.06 


6.34 


6.14±0.04 


6.49 


ICF(S+ +S 2+ ) 


1.12+0.10 


1.18 


1.50±0.15 


1.57 


1.61±0.13 


1.58 


12+log(S/H) 


6.27+0.12 


6.77 


6.49±0.11 


6.67 


6.48±0.09 


6.75 


12+log(Af+/H+) 


6.68+0.04 


6.62+0.03 


6.03 


5.94 


5.92±0.06 


5.98 


log(N/0) 


-1.06 ±0.10 


-1.12 ±0.08 


-1.11 


-1.16 


-1.26±0.12 


-1.32 


12+log(7Ve 2 +/i/+) 


7.27+0.03 


7.42+0.03 


7.30±0.06 


7.27±0.05 


7.33±0.03 


7.44±0.02 


log(Ne/0) 


-0.59+0.06 


-0.59±0.06 


-0.63±0.06 


-0.60±0.07 


-0.65±0.05 


-0.64±0.03 


12+log(Ar 2 +/H+) 


5.50+0.05 


5.67 


5.65±0.06 


5.62 


5.66±0.04 


5.70 


12+log(Ar 3 +/ J H'+) 


4.37+0.05 




4.32±0.07 


4.37±0.05 


4.36±0.05 


4.43±0.05 


ICF(Ar 2 ++Ar 3 +) 


1.21+0.06 




1.01±0.01 


1.02 


1.01±0.01 


1.02±0.01 


12+log(Ar/H) 


5.59+0.08 




5.66±0.06 


5.63 


5.66±0.05 


5.70 


12+log(Fe 2 +/fl'+) 


5.50+0.05 


5.66±0.07 


5.42±0.06 


5.46±0.06 


5.48±0.05 


5.59±0.05 


ICF(Fe 2 +) 


3.09+0.33 


3.52±0.32 


6.04±0.84 


6.63 


7.00±0.83 


6.72±0.63 


12+log(Fe/H) 


5.99+0.10 


6.20±0.11 


6.20±0.12 


6.28 


6.32±0.10 


6.42±0.09 


log(S 23 ) 


-0.22+0.02 




-0.02±0.03 




-0.10±0.02 






4600 4700 

Wavelength < A) 



Figure 6. Wolf-Rayct features in SDSS J162410.11-002202.5. 



have been derived from the values of the directly measured 
T([Sm]) and the near-IR [Sill] AA9069, 9532 A lines. The 
total sulphur abundance has been calculated using an ICF 
for S + ±S 2+ according to Barker's (1980 ) formula, which is 
based on the photoionization models bv lStasiriskal lll97ct) : 



ICF(S + + S 2 



N(Q+) 
N(O) 



1 -l/a 



where a = 2.5 gives the best fit to the scarce observational 
data on S 3+ abundances (Perez-Montero et al. 2006). 

The ionic abundance of nitrogen, N + /H + has been 
derived from the intensities of the AA 6548, 6584 A lines 
and the measured electron temperature of [Nil] for objects 
SDSS J002101.03+005248.1 and SDSS J162410. 11-002202. 5. 
Then, the N/O abundance has been derived under the as- 
sumption that 



JV 

o 



0+ 



In the case of SDSS J003218.60+150014.2, the [Nil] A 5755 A 
line could not be measured and therefore we have assumed 
that T([Nll]) is equal to T([Oll]) and we have derived di- 
rectly the N/O ra tio according to the expression given in 
IPagel et alj ll992l) . 

Neon is only visible in all the spectra via the [Nelll] 
emission line at A3868A. For this ion we have taken the 
electron temperature of [Oni], as representative of the high 
excitation zone. The total abundance of neon has been cal- 
culated assuming that: 



Ne 

~o 



Ne 2 + 

o 2 + 
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Izotov et al. (2004) point out that this assumption can lead 
to an overestimate of Ne/H in objects with low excitation, 
where the charge transfer between 2+ and H° becomes 
important. Nevertheless, in our case, it is probably justified 
given the high excitation of the observed objects. 

The main ionization states of Ar in ionized regions are 
Ar 2+ and Ar 3+ . The abundance of Ar 2+ has been calcu- 
lated by means of the emission line of [Arm] A 7136 A as- 
suming that T([Arm]) r*T([Siii]) (Garnett, 1992), while the 
ionic abundance of Ar 3+ has been calculated under the as- 
sumption that T([Anv]) « T([Om]) and using the emission 
line of [Anv] A 4740 A. We could have used the blended emis- 
sion line [Ariv]+Hel at A 4713 A subtracting the helium con- 
tribution. However, due to the relative abundance of these 
species and the signal-to-noise ratio of this blended line, we 
prefer not to estimate this ionic abundance with such a large 
error. The total abundance of Ar has been calculated using 
the ICF(Ar 2+ +Ar 3+ ) given by Izotov et al. (1994) which, in 
turn, has been derived from the photo-ionization models by 
Stasiriska (1990) as: 



ICF(Ar 2+ +Ar 3+ ) = 
O 



o+ 

'91 ( -q 



1.14 



O 



0.077 



91 

o 



+ \ 3-1 



Finally, for iron we have used the emission line of 
[Fein] A 4658 A and the electro n temperature of [Qui]. W e 
have taken the ICF(Fe 2+ ) from iRodri'guez fc Rubirl J2004I) . 
which yields: 



ICF(Fe 2 



(o 2+ ) 



1 



Q 2 + 

o+ 



The ionic and total abundances for each observed ele- 
ment are presented in columns 2, 4 and 6 of table [TU1 along 
with their corresponding errors. 



6 DISCUSSION 

6.1 Comparison with SDSS data 

In order to compare our results with those provided by SDSS 
spectra, we have measured the emission line intensities and 
equivalent widths on the SLOAN spectra of the three ob- 
served objects in the same way as described in Section 3. In 
order to allow an easy comparison between the results on 
both sets of spectra, we have listed these values in columns 
6 to 9 of tables H and El 

Strong emission line fluxes relative to H/3 measured 
on WHT and SDSS spectra differ by less than 10% for 
SDSS J003218.60+150014.2 and SDSS J162410. 11-002202. 5 
and about 30% for SDSS J002101.03+005248.1. This is 
partially compensated by differences in the derived red- 
dening constant so that reddening corrected emission line 
intensities relative to H/9 differ by less that 10% for 
the [On] AA 3727,29 A line, about 15% for the weak [Om] 
A 4363 A and [Sill] A 6312 A lines and only a few percent for 
the strong [Om] A 5007 A line. In fact, given the difference in 
aperture between both sets of observations, 0.5 and 3 arc- 
sec for WHT and SDSS respectively, some differences are to 
be expected. While we have probably observed the bright 
cores of the galaxies where most of the light and present 



star formation is concentrated, the SDSS observations map 
a more extensive area which could include external diffuse 
zones. This is evidenced by the more conspicuous underly- 
ing stellar population detected in the WHT spectra which 
leads to lower values of the emission line equivalent widths. 
The best agreement between the two sets of measurements 
is found for SDSS J003218.60+150014.2, which is probably 
the more compact object, as evidenced by the difference in 
the measured H/3 fluxes which is only a factor of 2 and the 
close agreement between the measured equivalent widths on 
the two spectra. 

The SDSS spectra have been analyzed following the 
same methodology as explained in Sections 4 and 5 for the 
derivation of temperatures and abundances although, due 
to the different nature of the observations: different spectral 
coverage, signal-to-noise ratio and spectral resolution, some 
further assumptions had to be made. These refer mainly to 
the temperatures of the different ions. Regarding sulphur, 
the SDSS data do not reach the 9000-9600 A range covered 
by the WHT spectra and therefore it was not possible to 
determine directly T([Sm]). The relation between T([Sm]) 
and T([Om]) is reproduced in Fig. |7| The sample used 
for comparison is a compilation of published data for which 
measurements of the nebular and auroral lines of [Om] and 
[Sill] exist, thus allowing the simultaneous determination of 
T([Om]) and T([Sm]). The sample is listed in Tables 1 and 
2 of Perez-Montero et al. (2006; objects not marked with 
a superscript b, as explained in the text), where the reader 
can find all the information concerning the data. The dashed 
line in the plot corresponds to the theoretical relation based 
on the grids of photoionization models described in Perez- 
Montero & Diaz (2005), 

r([5m]) = 1.05T([Om]) - 0.08 

which differs slightly from the semi-empirical relation by 
Garnett (1992) mostly due to the introduction of the new 
atomic coefficients for S 2+ from Tayal & Gupta (1999). The 
solid line in Fig. |7| corresponds to the actual fit to the data: 

T([Sm]) = (1.19 ± 0.08) T([Om]) - (0.32 ± 0.10) 

The individual errors have not been taken into account in 
performing the fit. This fit is different from that found 
by Garnett (1992) that seems to reproduce well the M101 
Hn region data analyzed by Kennicutt, Bresolin & Garnett 
(2003; KBG03). This is mostly due to the larger temperature 
baseline that we use. The object with the highest T([Sm]) in 
KBG03 is NGC5471A (12800 °K); our sample includes high 
excitation Hn galaxies with T([Sm]) up to 24000 °K, while 
including at the same time KBG03 sample. The introduc- 
tion of the high excitation objects make the relation steeper 
and increases the error of the calibration. This illustrates 
the danger of extrapolating relations found for a restricted 
range of values. We have used our empirical calibration in 
order to obtain T([Sm]) for the SDSS spectra. The estimated 
errors introduced by the calibration are of the order of 12% 
for TQSiii]), i.e. between 1400 and 1500 °K for the observed 
objects. 

Likewise, the SDSS spectrum of J003218. 60+150014.2 
does not include the lines of [On] at AA 3727,29 A and there- 
fore it was not possible to derive T([Oll]). Again, we have re- 
sorted to the model predicted relationship between T([On]) 
and T([Om]) found by Perez-Montero & Diaz (2003) that 
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2.6 - 



O HII galaxies 
2.4 - A GEHRs 




Tfloinj) 

Figure 7. Relation between T([Sm]) and T([Om]) for the ob- 
served objects (solid circles), and Hll galaxies (open circles), Gi- 
ant Extragalactic Hn regions (upward triangles) and diffuse Hll 
regions in the Galaxy and the Magellanic Clouds (squares) , for 
which data on the auroral and nebular lines of [Oin] and [So] 
exist (see Perez-Montero et al. 2006). The temperatures are in 
units of 10 4 °K. 



takes explicitly into account the dependence of T([0n]) on 
electron density. Three model sequences are represented in 
Fig.[S]correspondeing to three different values of the density: 
10, 100 and 500 cm -3 (dashed lines in the plot). The model 
sequence for n= 100 cm -3 is very similar to the one de- 
rived from the models by Stasiriska (1980). The sample used 
for comparison comprises the objects from Perez-Montero & 
Diaz (2005) for which the derivation of T([On]) and T([Om]) 
has been possible. In this case, due to the dependence of 
TQOll]) on electron density, there is not a single empirical 
calibration and it is not possible to give an estimate of the 
error introduced by the application of this procedure. 

In the case of nitrogen, the [Nil] A 5755 A line 
could be measured only in the SDSS spectrum of 
J002101. 03+005248.1, due to poor signal to noise. In the 
other two cases, the assumption T([Nn]) — T([Oll]) has been 
made. This assumption is usually made in standard analy- 
sis techniques; however, there are not enough data for Hll 
galaxies to test it empirically. 

Finally, the Balmer continuum temperature could not 
be calculated from the SDSS spectra due to lack of spectral 
coverage. 

Concerning abundances, the S 2+ /H + abundance ra- 
tios had to be calculated using the intensity of the weak 
auroral [Sin] line at A 6312 A and, in the case of SDSS 
J003218.60+150014.2, the 0+/H+ abundance ratio was de- 
rived using the [On] AA 7319,7330 A lines following the pro- 
cedure described by Kniazev et al. (2004). 

The values of electron density and temperatures, ionic 
and total abundances derived from the SDSS spectra are 
listed in columns 3, 5 and 7 of table 1101 along with their 
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Figure 8. Relation between T([On]) and T([Om]) for the ob- 
served objects (solid circles) and Hll galaxies (open circles), Gi- 
ant Extragalactic Hll regions (upward triangles) and diffuse Hll 
regions in the Galaxy and the Magellanic Clouds (squares) from 
Perez-Montero & Diaz (2005). The dashed lines correspond to 
photoionization models from Perez-Montero & Diaz (2003) for 
electron densities n e = 10, 100 and 500 cm -3 . The model se- 
quences from Stasiriska (1980; solid line) and Stasinska (1990; 
dashed-dotted line) are also shown. The temperatures are in units 
of 10 4 °K. 

corresponding errors in the cases where they have been de- 
rived from measured emission lines intensities. Otherwise, 
since the uncertainties introduced by the different assump- 
tions made and the theoretical models used are impossible to 
quantify, no formal errors are given. These quantities should 
be considered as estimates and be used with caution. 

For SDSS J003218. 60+150014.2 the values we have ob- 
tained for densities, temperatures and abundances from the 
WHT and SDSS spectra are in excellent agreement within 
the observational errors, as expected from the close agree- 
ment between the measured emission line intensities. For 
the other two objects, the agreement can be considered as 
satisfactory, taking into account the difference in aperture 
between both sets of observations. 



6.2 Gaseous physiscal conditions and element 
abundances 

For the three observed objects we have been able to mea- 
sure four electron temperatures: T([Om]), T([On]), T([Siii]) 
and T([Sn]). T([Nll]) has also been measured in two of the 
objects. The good quality of the data has allowed to reach 
accuracies of the order of 1% for T([Om]), 3% for T([On]) 
and 5% in the case of T([Nll]), T([Sll]) and T([Sm]). Figs. 03 
HOI and trTl show the range of the measured line temperatures 
and electron density for objects SDSS J002101.03+005248.1, 
SDSS J003218.60+150014.2 and SDSS J162410. 11-002202. 5, 
respectively. The data for the latter two are very simi- 
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lar, although T([Nn]) could not be measured for SDSS 
J003218.60+150014.2. The width of the bands correspond 
to one <r error. In the two cases T([Oiii]) and T([Sm]) over- 
lap to some extent, while T([0n]) is lower than T([Oiii]) in 
the first case and slightly higher than T([Om]) in the second. 

The [On] and [Oiii] temperatures of the observed ob- 
jects are shown in Fig. |H| together with the data of Hn 
galaxies, and galactic and extragalactic Hn regions from 
Perez-Montero & Diaz (2005) and the photoionization mod- 
els from Perez-Montero & Diaz (2003) for electron densities 
n e = 10, 100 and 500 cm -3 . The model sequences from 
Stasfiska (1980; solid line) and Stasihska (1990; dashed- 
dotted line) are also plotted. This figure shows the effect 
of the electron density on the derived [On] temperature 
which decreases with increasing density. The data points 
populate the region of the diagram spanned by model se- 
quences but the observational errors are too large to al- 
low an adequate test of the photoionization models them- 
selves. The temperatures derived for two of the observed 
objects: SDSS J003218.60+150014.2 and SDSS J162410.ll- 
002202.5 lie close to the theoretical relation for an electron 
density of n e = 10 cm _3 , while for the other one: SDSS 
J002101. 03+005248.1, they lie closer to sequences for elec- 
tron densities larger than 100 cm -3 . This is consistent with 
the trend shown by the values of n e measured for the three 
objects. 

For this last object, the value of T([On]) derived from 
the measured T([Oiii]) according to the widely used fit to 
the photoionization models by Stasinska (1990), is larger 
than the measured one by 2200 °K which translates into a 
lower + /H + ionic ratio by a factor of 4 and a lower total 
oxygen abundance by 0.17 dex. It should be remarked that 
this is the procedure usually followed, even in the cases in 
which measured values of T([Oll]) exist (see Izotov, Thuan 
& Lipovetsky 1994). Also, in doing so, no uncertainties are 
attached to the T([On]) vs. T([Oni]) relation with the final 
outcome of a reported T([On]) which carries only the usually 
small observational error found in the derivation of T([Oiii]) 
and translates into very small errors in the oxygen ionic and 
total abundances. Thus it is possible, and even frequent, to 
find reported values of T([On]) and T([Sm]) with quoted 
fractional errors lower than 1% and absolute errors actuall y 
less than that quoted for T([Om]) Jlzotov fe Thuanlll99ct . 
and ionic + /H + ratios with errors of only 0.02 dex. 

Recently, thi s procedure has been justified by 
llzotov et alj (120061 ') on the basis of an analysis of SDSS data. 
In this work it is argued that, despite a large scatter, the 
relation between T([0n]) and T([Om]) derived from obser- 
vations follows generally the one obtained by models, and 
therefore they adopt T([Oll]) as derived from their photoion- 
ization models. However, the large errors attached to the 
electron temperature determinations, in many cases around 
± 2000 °K for T([On]), actually precludes the test of such a 
statement. In fact, most of the data with the smallest error 
bars lie below and above the theoretical relation. 

Regarding the [Sill] temperatures, our values fit well the 
relation between T([Oiii]) and T([Sm]) found by Garnett 
(1992) as can be seen in Fig. [7] For the objects for which we 
have measured T([Nll]), this is equal to T([Sm]) within the 
errors. 

The abundances derived for the observed objects 
show the characteristic low values found in strong 
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Figure 9. Measured line temperatures and electron density for 
object SDSS J002101. 03+005248.1. The width of the bands cor- 
respond to one cr error. 
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Figure 10. Measured line temperatures and electron density for 
object SDSS J003218.60+150014.2. The width of the bands cor- 
respond to one cr error. 



line Hn galaxies (iTerlevich et al] Il99ll : iHovos fc Di'azl 

12006ft : 12+log(0/H) =8.0, within the errors, for SDSS 
J002101. 03+005248.1 and SDSS J003218.60+150014.2 and 
slightly lower 12+log(0/H) = 7.9 for SDSS J162410.ll- 
002202.5. The three objects have previous oxygen abun- 
dance determinations. They are part of the first edition 
of the SDSS Hn galaxies w i th oxy gen abundance catalog, 
presented by iKniazev et al] i2004ft . These authors derived 
total oxygen abundances of 12+log(0/H) =8.18+0.04 for 
SDSS J002101.03+005248.1, 12+log(0/H) = 8.07+0.02 for 
SDSS J003218.60+150014.2 and 12+log(0/H = 8.17+0.01 
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Figure 11. Measured line temperatures and electron density for 
object SDSS J162410.11-002202.5. The width of the bands corre- 
spond to one cr error. 



for SDSS J162410.11-002202.5, higher than ours by 
0.17 dex, but consistent with the values we obtain from 
the analysis of the SDSS spectra for two of the objects: 
SDSS J002101.03+005248.1 (12+log(0/H) = 8.19±0.03) 
and SDSS J162410. 11-002202.5 (12+log(0/H = 8.14±0.02). 
For SDSS J003218.60+150014.2, our analysis of its 
SDSS spectrum yields a total oxygen abundance 
12+log(0/H) =7.93±0.03, lower t han theirs by . 14 dex 
and closer to the value derived bv lUgrvumov et alJ (l2003f) 
in their Hamburg / SAO Survey. The spectrum analyzed 
by iKniazev et"aiT J2004I) was extracted from the first 
SDSS data release. This might point to a difference in the 
calibration routines of SDSS spectra from one release to the 
other. 

The logarithmic N/O ratios found for the two galaxies 
for which there are T([Oll]) and T([Nll]) determinations are 
-1.06±0.10 and -1.26±0.12. For the third object we estimate 
a log(N/0) ratio of -1.11. They point to a constant value 
within the errors. It is worth noting that an analysis of the 
data along the lines discussed above, i. e. T([Oll]) derivation 
from T([Om]) according to Stasihska (1990) models, would 
provide for one of the objects, SDSS J002101.03+005248.1, 
an N/O ratio larger by a factor of 3 (log N/O =-0.64). 
Therefore it is not unplausible that part of the scatter found 
in the N/O vs. O/H diagram may be due to the methodol- 
ogy employed in the derivation of the N/O ratios. This is an 
important effect that should be explored further. 

Finally, the log(S/0) ratios found for the three 
objects are: -1.83±0.16, -1.44±0.14 and -I.57i0.ll, 
barely consistent with solar (S/0 = -1.39) except for 
J002101. 03+005248.1 which is lower by a factor of about 
2.7. The analysis of the data according to the conventional 
lines described above would rise this S/O ratio somewhat, 
but still keeping it below the solar value. 

On the other hand, the measurement of the [Sin] IR 



lines allows the calculation of the su lphur abundance pa- 
rameter S 23 (IVilchez fc Estebadll996h . 

_ [511] AA 6717. 31 + [Sin] AA 9069, 9532 

523 w 

This parameter constitutes probably the best empir- 
ical abundance indicator for Hll galaxies, since contrary 
to what happens for t he widely used O 23 (R23) parame- 
ter iPaeel et all Il979l: lAUoin et alJll979l) . the calibration 
is linear up to solar abundances, thus solving the degen- 
eration problem usually presented by this kind of objects. 
This is particularly dramatic for objects with log023 > 0.8 
and 12+log(0/H) >8.0 which can show the same value of 
the abundance parameter while having oxygen abundances 
that differ by up to an order of magnitude. About 40% 
of the observed Hll galaxie s belong to this category (see 
iDfaz fc Perez-Monterolll99flt) . The logarithm of the S23 pa- 
rameter derived for our objects are given in the last row 
of table 1101 Fig. 1121 shows the points corresponding to the 
objects in the logS23 vs. 12+log(0/H) diagram, together 
with their observational error, and labeled as WHT objects. 
The rest of the data correspond to Hn galaxies from Perez- 
Montero & Diaz (2005) tog ether with the data of ten SDSS 
BCD gal axies presented bv | K niazev e t alJ (120041) and ana- 
lyzed bv IPerez-Montero et alJl200^W these objects are la- 
beled as SDSS objects). The solid line shows the calibration 
by Perez- Montero & Diaz (2005). Three objects are seen 
to clearly deviate fro m this line. Two of th em correspond to 
two of the galaxies of lKniazev et al.l <2004h for which no data 
on the [On] 3727 line exist and therefore the + /H + ratio 
is derived from the red [On] 7325 lines. Both objects show 
low ionization parameters as estimated from the [Sn] / [Sin] 
ratio and relatively high values of O/H as derived from 
the [Nil] /Ha calibration (|Denicol6. Terlevich fc Terlevichl 
12002ft . The third object is Mrk709, an object that shows 
similar characteristics (see Perez-Montero & Diaz 2003). 
These objects might be affected by shocks and deserve fur- 
ther study. Actually the accuracy of the S23 calibration for 
Hll galaxies, as a family, is only 0.10 dex (see Perez-Montero 
& Diaz 2005). But obviously, more observations are needed 
in order to improve the S23 calibration and truly understand 
the origin of the observed dispersion. 

6.3 Ionization structure 

The ionization structure of a nebula depends essentially on 
the shape of the ionizing continuum and the nebular ge- 
ometry and can be traced by the ratio of successive stages 
of ionization of the different elements. With our data it is 
possible to use the + /0 2+ and t he S + /S 2+ to probe the 
nebular ionization structure. In fact. lVilchez fc Paeell (Il988l) 
showed that the quotient of these two quantities that they 
called "softness parameter" and denoted by rj is intrinsi- 
cally related to the shape of the ionizing continuum and 
depends on geometry only slightly. If the simplifying as- 
sumption of spherical geometry and constant filling factor 
is made, the geometrical effect can be represented by the 
ionization parameter which, in turn can be estimated from 
the [Oii]/[Oiii] ratio. Actually, the [Oii]/[Om] ratio depends 
on stellar effective temperature which, in turn, depends on 
metallicity, in the sense that for a given stellar mass, stars of 
higher metallicity have a lower effective temperature. Since 
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Figure 12. Relationship between logS23 and metallicity, repre- 
sented by 12+log(0/H), for the observed objects (solid circles), 
the Hn galaxies from Perez-Montero & Diaz (2005) with data on 
the [Sin] AA 9069,9532 A emission lines (open circles), and SDSS 
BCDs galaxies (asterisks) from Perez-Montero et al. 2006. The 
solid line shows the calibration made by Perez-Montero & Diaz 
(2005). 



the three observed objects have similar values of this quo- 
tient (between 0.2 and 0.3), and show oxygen abundances 
in a very narrow range, we can assume they also share a 
common value for their ionization parameter. Under these 
circumstances, the value of rj points to the temperature of 
the ionizing radiation. 

In panel (a) of Fig. M 

we show the relation between 
log(0+/0 2+ ) and log(S+/S 2+ ) for the observed objects 
(solid circles for values calculated using our methodology 
and solid diamonds for values derived using the conventional 
method) and the Hn galaxies from Perez-Montero & Diaz 
(2005) for which the different oxygen and sulphur ionic ra- 
tios can be derived (open circles). In this diagram, diago- 
nal lines correspond to constant values of 77; the solid line 
corresponds to log 77 = 0.00. The region above this line cor- 
responds to 77 > 1 and the region below the line corresponds 
to r) < 1. The dashed lines corresponds to log 77 = -0.35 and 
log 77 = 0.2 . As can be seen in the figure, Hn galaxies occupy 
the region of the diagram with -0.35 < log 77 < 0.20, which 
corresponds to high values of the ionizing temperature ac- 
cording to Vilchez & Pagel (1988). Two of the observed ob- 
jects lie on the log 77 = -0.35 line. The third one, the one with 
the highest excitation, has log 77 = -0.15. An analysis along 
the conventional lines commented in the previous subsec- 
tion, would yield the same value of 77 for the highest excita- 
tion object (although with different values of both + /0 2+ 
and S + /S 2+ ratios) but widely different values of 77 for the 
other two objects: 0.19 and 0.30 indicating much lower ioniz- 
ing temperatures. This is inconsistent, however, with what 
is found if the same exercise is performed using the cor- 
responding quotients of emission lines: log([On]/[Oni]) vs 



log([Sll]/[Sni]), which does not require explicit knowledge 
of the line temperatures involved in the derivation of the 
ionic ratios, and therefore does not depend on the method 
to derive or estimate these temperatures. Panel (b) in Fig. 
H3l shows then the purely observational counterpart of panel 
(a). In this diagram diagonal lines represent constant values 
of log 77' defined by Vilchez & Pagel (1988) as: 

[On] A A 3727, 29 / [Om] A A 4959, 5007 



[Sn] AA 6717, 31 / [Sin] AA 9069, 9532 
0.14 



0.16 



where t e is the electron temperature in units of 10 . 77 
and 77' are related through the electron temperature but 
very weakly, so that a change in temperature from 7000 
to 14000 °K implies a change in log 77 by 0.1 dex, inside ob- 
servational errors. Also, log 77' is always less than log 77. In 
this diagram, the data for the three objects indicate values 
of log 77' between -0.39 and -0.25, consistent with what is 
found by our analysis but inconsistent with the results ob- 
tained following the conventional methodology. Metallicity 
calibrations based on abundances derived according to this 
conventional method are probably bound to provide metal- 
licities which are systematically too high and should there- 
fore be revised. 



6.4 The temperature fluctuation scheme. 

At th e end of the 60s and beginn ing o f the 70s, IPei mbertl 
Jl967fl . |Peimbert fc Costercl (Il969l) . and lPeimbertJ Jl97ll) es- 

tablished a complete analytical formulation to study the dis- 
crepancies between the abundances relative to hydrogen de- 
rived from recombination lines (RLs) and from collisionally 
excited lines (CELs) when a constant electron temperature 
is assumed. In the first of these works, Peimbert proposed 
that this discrepancy is due to spatial temperature varia- 
tions which can be characterized by two parameters: the 
average temperature weighted by the square of the density 
over the volume considered, To, and the root mean square 
temperature fluctuation, t 2 . They are given by 



T pr+) = 



and 



t 2 {X i+ ) 



[T e N e N(X t+ )dV 
f NeNiX^dV 

J(T e - T (X i+ )f N e N(X i+ ) dV 
ToiX^) 2 J N e AT(X*+) dV 



(1) 



(2) 



where N e and N(X +l ) are the local electron and ion den- 
sities of the observed emission lines, respectively; T e is the 
local electron te mperature; and V is the observed volume 
JPeimbertlll967l) . 

It is possible to obtain the values of To and t 2 using dif- 
ferent methods. One possibility is to compare the electron 
temperatures obtained using two independent ways. Gener- 
ally, temperatures originated in different zones of the nebula 
are used, one that represents the warmest regions and an- 
other representative of the coldest ones. In the absence of 
temperatures derived from recombination lines, the temper- 
atures estimated from the hydrogen discontinuities, either 
from Balmer or Paschen series, representative of the tem- 
perature of the neutral gas, can be used. 
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log (S + /S 2+ ) log ([SII]U6717,31/[SIII]U9069,9532) 

Figure 13. Panel (a): Relation between log(0 + /0 2 ~'~) and log(S~'~/S 2 ~'~) for the observed objects (solid circles for values calculated 
using our methodology and solid diamonds for values derived using the conventional method) and the Hll galaxies from Perez-Montero 
&; Diaz (2005) for which the different oxygen and sulphur ionic ratios can be derived (open circles). Diagonals in this diagram correspond 
to constant values of r\. Panel (b): Relation between log([Oll]/[Om]) and log([Sn]/[Sin]) for the observed objects (solid circles) and Hll 
galaxies from Perez-Montero & Diaz (2005) with data on the [Sill] AA 9069,9532 A emission lines (open circles). Diagonals in this diagram 
correspond to constant values of r\\ 



We have followe d iPeimbert et"aH feOOd 120021. 12004 
and lRuiz et alJ i200otl to derive the values of To and t 2 for 
our WHT spectra by combining the Balmer temperature. 
T(Bac), and the temperature derived from the collisional 
[Om] lines, T([Om]), therefore assuming a simple one-zone 
ionization scheme. Then, the relation for these two temper- 
atures is given by: 



T([Om]) = To 
and 



1 / 91300 
+ 2 I To 



3)f 



T(Bac) = T 1 - 1.67* 



(3) 



(4) 



The solution of this system of equations for each galaxy 
along with their corresponding errors are listed in table 
1111 The temperature fluctuations are almost negligible for 
two of our objec ts. They have t 2 - values very similar to the 
ones derived by iLuridiana et al.l i2003l) by combining ob- 
servations from the literature and photoioni zation models 
for som e BCDs and extragalactic Hn regions. iGuseva et alJ 
i2006l) used the Balmer and Paschen jumps to determine 
the temperatures of the H + zones of 22 low-metallicity Hll 
regions in 18 BCD galaxies, one extragalactic Hn region 
in M101 and 24 Hll emission-line galaxies selected from 
the DR3 of the SDSS. They found that these temperatures 
do not differ, in a statistical sense, from the temperatures 
of the [Oni] zones, given t 2 -values are close to zero. The 
greater t 2 -value obtained for SDSS J003218.60+150014.2 is 
in the range of the values derived for giant e xtragalactic Hll 
regions (see iGo nzalez- Delgado et aljTl994l : IPeimbert et alJ 
120051: Ijamet et alj|2005l) or Hll regions in the Magellanic 



Table 11. To and t parameters 



name 


To 


t 2 


SDSS J002101. 03+005248.1 


1.24±0.35 


004+ ' 044 


SDSS J003218.60+150014.2 


1.08±0.21 


0.066±0.026 


SDSS J162410.11-002202.5 


1.24±0.30 


001 + 037 
u.uiu_ 001 



To in 10 °K. Note that t is always greater than zero. 



Clouds, 



such as 3 
aH l2003h . 



Doradus, LMCN11B and SMCN66 



In the limit of low densities and small optical depths, 
and for t 2 much lower than one, the electronic temperature 
for helium, T(Hell), is proportional to (a) and /3, the av- 
erage value of the power of the temperature for the helium 
lines used to calculate the ionic abundances of He + , and the 
corresponding one for H/3, respectively. Then, for (a) differ- 
ent from f3, this temperature is given by equation (14) of 
Peimbert (1967): 



T(Hen) = T{Heii,Hu) = T 



1+ ((«>+/?-!)- 



(5) 



The value of the power of the temperature for each 
heliu m line in the low den sity limit has been obtained 
from iBeniamin et al.l il999T) . We have calculated (a) as 
the average value of a weighted accor ding to the ob- 
servational errors iPeimbert et alJ |2000D . We have ob- 
tained values of (a) equal to -1.37, -1.42, and -1.43 for 
SDSS J002101.03+005248.1, SDSS J003218.60+150014.2, 
and SDSS J162410. 11-002 202. 5, respectively The v alue of 
P has been obtained from IStorev fc Hummeil dl99fJ) and is 
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equal to -0.89. The results for T(Hell) and their correspond- 
ing errors are listed in the last row of table [S] 

The temper ature for H/3, T(H/3), c an be calculated from 
equation (20) in lPeimbert fc Costerol (^969) as: 



T{H(3) = To 



1-0.95 V 



(6) 



where we have taken (3 = -0.89 as above. The derived values 
for T(H/3) and their errors are also given in table |S] 

The line temperature for a collisionally excited line, 
CEL, for t 2 < 1, ( AEcEL/kT -l/2) ^0 and a / is given 
by equation (20) of lPeimbert fc Costerol (ll96ST) : 



Tci 



Toil 



+ 



(AEcEL/kTof - 3AE C EL/kTa + 3/4 



AEr 



-JkT - 1/2 



(7) 



where AEcel = A_B m „ is the energy difference, in eV, be- 
tween the upper (m) and lower (n) levels respectively of the 
atomic transition that produces the line. 

For the case of t 2 > 0, and assuming a one-zone ioniza- 
tion scheme, we can derive the ionic abundances using the 
values calculated for t 2 equal to zero, and equation (15) of 
iPeimbert et~aH J20041 : 



N(X +1 ) 



x exp 



T{Ht3)-°* 9 T{\ mn f 
T([Om])-°- 37 

AE ae 

+ 



fcTQOnrj) kT(\ mn ) 



7V(X H 



(8) 



where T(A mn ) is given by equation Using this expres- 
sion, we have calculated the effect of the temperature fluc- 
tuations on the 2+ /H + abundance. In fact, given the high 
excitation of the object, this ionic abundance carries the 
highest weight in the total abundance of oxygen. The recal- 
culated value is 12+log(0 2+ /H + ) = 8.09±0.12 which yields 
a total oxygen abundance 12+log(0/H) = 8.12±0.11. These 
values are higher than those given in Table [TU1 by 0.23 and 
0.22 dex respectively. 



7 CONCLUSIONS 

We have performed a detailed analysis of newly obtained 
spectra of three Hn galaxies selected from the Sloan Digital 
Sky Survey Data Release 2, covering from 3200 to 10550 A 
in wavelength. For the three objects we have measured four 
line temperatures: T([Om]), T([Sm]), T([On]) and T([Sn]) 
and the Balmer continuum temperature T(Bac). For two of 
the objects we have also measured T([Nll]). These measure- 
ments and a careful and realistic treatment of the observa- 
tional errors yield total oxygen abundances with accuracies 
between 5 and 9%. The fractional error is as low as 2% 
for the ionic 2+ /H + ratio due to the small errors associ- 
ated with the measurement of the strong nebular lines of 
[Om] and the derived T([Oiii]), but increases to 15% for the 
0+/H+ ratio. The accuracies are lower in the case of the 



abundances of sulphur (of the order of 20% for S + and 12% 
for S 2+ ) due to the presence of larger observational errors 
both in the measured line fluxes and the derived electron 
temperatures. The error increases further (up to 30%) for 
the total abundance of sulphur due to the uncertainties in 
the ICF. 

This is in contrast with the small errors quoted for line 
temperatures other than T([Om]) in the literature, in part 
due to the commonly assumed methodology of deriving them 
from the measured T([Oni]) through theoretical relations. 
These relations are found from photoionization models and 
no uncertainty is attached to them; therefore, the obtained 
line temperatures carries only the observational error found 
for the T([Oni]) measurement. If this methodology were to 
be adopted, the theoretical relations should be adequately 
tested and empirical relations between the relevant line tem- 
peratures should be obtained in order to quantify the corre- 
sponding model uncertainties. 

We have compared our obtained spectra with those 
downloaded from the SDSS DR3 finding a satisfactory agree- 
ment. The analysis of these spectra yields values of line tem- 
peratures and elemental ionic and total abundances which 
are in general agreement with those derived from the WHT 
spectra, although for most quantities, they can only be taken 
as estimates since, due to the lack of direct measurements of 
the required lines, theoretical models had to be used whose 
uncertainties are impossible to quantify. Unfortunately, the 
spectral coverage of SDSS precludes the simultaneous ob- 
servation of the [On] AA 3727,29 A and [Sm] AA 9069, 9532 A 
lines, and therefore the analysis can never be complete. 

The ionization structure found for the observed objects 
from the + /0 2+ and S + /S 2+ ratios points to high val- 
ues of the ionizing radiation as traced by the values of the 
"softness parameter" 77 which is less than one for the three 
objects. Line temperatures derived from T([Oni]) with the 
use of current photoionization models yield for the two high- 
est excitation objects, much higher values of r\ which would 
imply lower ionizing temperatures. This is however inconsis- 
tent with the ionization structure as probed by the measured 
emission line intensities. 

Finally, we have measured the Balmer continuum tem- 
perature for the three observed obje cts and der i ved th e tem- 
perature fluctuations as defined bv iPeimbertl il967f) . Only 
for one of the objects, the temperature fluctuation is signif- 
icant and could lead to higher oxygen abundances by about 
0.20 dex. 
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